Wigner distributions for finite state systems without redundant phase point operators 



o 
o 

(N 



S. ChaturvedE 
School of Physics, University of Hyderabad, Hyderabad 50004-6 

N. Mukundifl 

Centre for High Energy Physics, Indian Institute of Science, Bangalore 560012 

R. SimoiH 

The Institute of Mathematical Sciences, C. I. T. Campus, Chennai 600113 

We set up Wigner distributions for A'' state quantum systems following a Dirac inspired approach. 
In contrast to much of the work on this case, requiring a 2A'' x 2A'' phase space, particularly when 
A'^ is even, our approach is uniformly based on an A'' x A'^ phase space grid and thereby avoids 
the necessity of having to invoke a 'quadrupled' phase space and hence the attendant redundance. 
Both A'^ odd and even cases are analysed in detail and it is found that there are striking differences 
between the two. While the A^ odd case permits full implementation of the marginals property, the 
^1^' even case does so only in a restricted sense. This has the consequence that in the even case one is 

D , led to several equally good definitions of the Wigner distributions as opposed to the odd case where 

the choice turns out to be unique. 
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^ ■ I. INTRODUCTION 

The Wigner distribution, as originally defined by Wigner [3| for the case of a Cartesian configuration space has, 
over the years, been extended and generalised in many directions [3|-[23j. These extensions and generalisations (most 
K*" ' of which have been elegantly unified in a recent work of Ferrie and Emerson [l^l using the language of frame theory) 
include situations where the coordinates take values in 
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m ■ .a finite field F„ 



• a, ring 
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a finite group (abelian or non abclian) 
• the manifold of a semi simple compact Lie group 

" The first two are of particular interest owing to potential applications to quantum information processing and quantum 
state estimation. In the present work we focus our attention on the second case and set up the phase space and Wigner 
^ . distributions thereon using what we call a Dirac inspired approach to Wigner distributions [20| . This approach 

■ ■ ■ ' essentially requires computing the square root of a certain kernel which brings with it undetermined signs, one at 
each phase point. We examine the question to what extent these signs can be fixed or related to each other by the 
marginals property i.e by demanding that the phase space averages of the Wigner distribution along an 'isotropic 
line' and the lines 'parallel' to it yield a probability distribution. For the case when N is odd one finds that one 
can consistently impose marginals property on all 'isotropic lines' which in turn uniquely fixes all the signs. The 
Wigner distribution thus obtained turns out to be the same as that already known in the literature. For the even 
case, on the other hand, this is not so. We find that the marginals property can not be consistently imposed on all 
the 'isotropic lines' but only on specific subsets thereof (i.e. on individual orbits under S'L(2,Zjv)). Therefore the 
best one can do is to demand the marginals property on the largest such subset. This is actually good enough as the 
'isotropic lines' in this subset ( and only in this subset) cover all phase points. We explicitly carry this out and find 
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that, unlike the odd case, not all signs get fixed and therefore one has a family of Wigner distributions characterised 
by different choices for the signs , all of which respect the restricted marginals property. We emphasise that here we 
work all along with a,n N x N phase space lattice and the associated phase point operators. This is in contrast 
to other formalisms to be found in the literature 0, Q where a 2N x 2N phase space lattice is invoked to arrive 
at a satisfactory definition of Wigner distributions for an iV-state quantum system. An interesting recent work by 
Bar-on [2^ based on symmetric informationally complete projection operator valued measures (SIC-POVMs) [l^l and 
a block design theory [25] inspired picture of phase space, reduces the number phase points from 47V^ to N'^ + N, but 
still has N phase points more than those required in the present work. 

A brief outline of this work is as follows: In Section II we set up our notation and summarise the properties of the 
operators used later. In Section III we briefly recapitulate the Dirac inspired approach to Wigner distributions and 
give the expressions for the phase points operators in a convenient form along with the conditions on the undetermined 
signs that appear therein in order that the phase point operators satisfy the standard marginals property. In Section 
IV we briefly summarise the properties of the isotropic lines in our phase space as made available in a recent work by 
Albouyjl^. In Section V we discuss the odd case and in Section VI summarise our result for the case when is a 
prime power. The case of general A^ is discussed inSection VII. In Section VIII we discuss the question of tomographic 
reconstruction for the case when N = 2^. We conclude with Section IX where we note that although for a given even 
N there is a great variety of Wigner distributions characterised by different choices for the unfixed signs, the actual 
number of choices available, in so far as the eigenvalues of the relevant phase point operators is concerned, is rather 
small. 



II. PRELIMINARIES 



Consider a quantum system described by a complex Hilbert space Ti. of dimension A^. Denote by {!<?)}, <1,P S 

Zjv two orthonormal bases - coordinate and momentum bases - related to each other by a finite Fourier transform : 

{q'\q) - <5g',«; (p'Ip) = (qIp) - ^"VVa; ^ = e^-/^ (i) 

Here Z^v stands for the ring of integers {0, 1, 2, • • • , A — 1} with addition and multiplication modulo A^. We use the 
notation [n] to denote n modulo A. On the Hilbert space H we introduce the familiar Weyl operators: 

C/ = ^4= ^ u;1\q){ql U \p) = \[p + 1]) ; U^U = U^ = l. 
q€ZN 

C/p = c.p« = ^ \q){q\ = ^ E (2) 

v = u;-p^Yl ^'"IpM^ V\q)^\[q + l]y, yty^yiV^i 
peZN 

V^^iu-^P^ J2 u;-P'^\p)ip\; |P)(P|-^ E 

PGZn qeZfi 

C/Py« = LuP^vmP = r'^PivmP, (3) 

where r = e^^l^ . We denote by Fq the discrete 'classical' phase space of A^ points a — {q,p) equipped with the 
symplectic product {cr,(j') — pq' — qp' between any two phase points a and a' . We denote by K, the Hilbert space of 
square summable functions on Fq: 

JC — i^(Fo) = Af^— dimensional Hilbert space 

= {/(a)eC|aero, || / f = E 1/(^)1'}- (4) 



Fq can be viewed as an abelian group of order A^ 
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(7 e Fo — > group element g{(j) = a'^W , 
generators a,h : ah — ha, a^ = h^ ~ e 

5((0,0)) = identity e; (5) 
composition : g{cj')g{(T) = g{{[q + q], [p + p])) = g{[(T + cr']); 

inversion: g^ay^ ^ g{{\N - q],[N - p])) = g{\N - a]). (6) 
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Thus Tq is the direct product Go x Gq where Gq is the (abehan) cyclic group of order N. 

Tq has N'^ distinct incquivalent one-dimensional unitary irreducible representations (UIR's), with characters la- 
belled by points cr' G Fq : 

a'eTo: UIR a ^ uj-p' , b ^ uji' ; 
Orthogonality : xl"X<y' = X*a"{(^)Xa'{<y) 

= tjP"Q-<l"P+l'p-p'q ^ N^Sa>>,„' (7) 

The set of functions {jfXa'i'^)} on To forms an orthonormal basis (ONB) for K,. The trivial UIR is a' = 0, Xo(o') = 1; 
and hence 

creTo aero 

Finally we introduce the displacement operators D{a) on 7Y, one for each phase point, and summarise their properties: 

(tgFo: D{a) = D{q,p)=T''PVmP = t-^pUPV''; (9) 
Unitarity : D{a)^D{a) = 1 onH. (10) 
Inverses : D{a)-^ = 7]aD{[N - a]), (11) 

1 for (7 = or p = 
(_1)9+P+W for 1 < < AT - 1 

^ ^-lN-q]lN-p]+qp ^^2) 

Composition : D{<t')D{(j) ^ t^"' ^"^{(1' , a)D{[(7' + a]), (13) 

1 for ct' + (T e Fo 

(-l)«'+« for q' + g < TV - l,p' +p > iV 
(_1)P'+P ioY q' + q> N ,p' + p < N - 1 
(^^ly' +q+p' +P+N for q' + q> N,p' +p> N 

_ ^{q +q){p' +p)-lq' +q]lp' +p] ^^4) 

Trace orthogonality : Tr{D{(T')^ D{<t)) ^ NS^' ,a- (15) 

The set of operators {^^D{a)} constitutes a complete irreducible trace orthonormal set of operators on H satisfying 
the relations: 

D{a')D{a)^Lu''^''^^D{a)D{a') (16) 
D{a')D{a)D{a')-^ ^Lo''^''''^D{a) (17) 

Thus {D(a)} form an irreducible unitary N dimensional ray representation of Fo on Ti.. Some useful relations are 
given below: 

a e Fo : T-^PDia) = Vm^ = ^ ^^"^'^'^ \p')ip'W){<l'l (18) 

\p)ip\ci){q\ = ]^ E '^^"^^'^r-^'^'Di'^')- (19) 
cr'ero 



e(cr', (t) = e(fT, cr') 



III. DIRAC INSPIRED APPROACH TO WIGNER DISTRIBUTIONS 

The central idea in this approach, discussed in detail in [20], is to initially associate with any operator A a phase 
space function A{a) constructed out of its mixed matrix elements such as (g|A|p) in such a way that the trace of the 
product of two operators is expressed as a phase space average of a kernel times the product of their phase space 



4 



functions, and then try to transform away this kernel. As shown in [20'|, this exercise entails finding a symmetric 
square root of the kernel K{a,a') defined below: 

K{a;a') ^ K{a';a) ^ uj^'^-'^">^P-P'l (20) 

A. The Kernel K and its symmetric square roots ^ 

The kernel K defines an operator on JC. The characters {xa' (c)} of To form a complete orthogonal set of eigenvectors 
of K in /C: 

J2 K{<T- o')xa" {a') = iVu;«"f "x." (a) (21) 

Using the results above for the eigenvalues and eigenvectors of K, a general square root ^ of ii" is defined by 

^X.„ ^ VNT^"p"s{a")x.", S{a") = ±1, a" e To. (22) 
At this point, there are N'^ sign choices. This freedom will be reduced as we proceed. The kernel ^(cr; a') is 

cr"ero 

= 1^ E r''"p"Sia")J^^^"^~^^''^"^ (23) 
<y"ero 

K{a; a') is symmetric under a a' . Demanding the same for ^(cr; a') places conditions on the signs 5((t"). Since in 
general 

E E fiW-l]) etc, (24) 



c"ero 



'?"p".';'('rr"V,,('^.'^")-('^'.'^"> 

((T ,(T )-(o-,(T ) 



o-"ero 

= ^ ^[JV-9"][JV-p"]5,(j^_^,,])^(a,<r">-(a'.<r"> 

o-"ero 

^ ^(ff) = r[^"«l[^^Pl-«P5([iV - a]) = ?7,5([A^ - a]). (25) 
So the number of independent 5((t)'s is about halved. We hereafter assume ^(cr;cr') is symmetric. 

B. Phase point operators W{a) 

For any choice of the square root kernel ^(cr; cr'), we define 

aero: W{a) ^ Vn Yl ^('t; a')b')(p'k')(9'l 

= E ^^"'^'^5(a')i^(a') (26) 

The condition (|25|) on S{(7) arising from the symmetry of ^ ensures hermiticity of W{a): 

^(^)^ = ^ E ^^^(^')i^(^')-' 



1 5] c.<-'-'^)^([A^-cT'])i?([iV-a']) 
(T'ero 

W{<7) (27) 
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Trace orthogonality of D{a) leads to 

Ti{W{a')W{a)) = NSa',a (28) 

Further Tr(L>(CT)) = N6a,o gives 

Tr(t?(a)) = 1, (29) 

So both {-^W{a)} and {:^-^(''')} ^'''6 trace orthonormal complete sets of operators on H. 
From the conjugation relations (|17p for D's we get: 



i?(a')l^^((T)I?(a')"' = ^ ^ '^^"^""^Sia")Lu''^''^"^Dia") 

<T"ero 

= M?([a + a']) (30) 

Recovery of standard marginals fixes some S{a): 

1 J2 = ^2 E c^^''-'^'5(g',p')^('z',p') 



peZjv q',p',p€ZN 



TV 



= |9)(g|^5(0,p') = l; (31) 
1 ^ t?(«,p) = \p)ip\ ^ Siq\0) = 1 (32) 



These are consistent with ([25|) . So at this stage 

5(g,0) = 5(0,p) = l; ^(a) = r;,5([7V-a]) (33) 
Some useful formulae are given below: 

Tr{W{a)D{a')^) = lu'^^''''^ S{a'), 
Ti{W{a)D{a')) = lu<^' '''^ S{a'), 

S{a)D{a) = ^ E ^^"^^'^Wia') (34) 

IV. ISOTROPIC LINES AND FURTHER MARGINALS 

The conditions so far on S{a) are given above in p3p . To generate more conditions, we consider more marginals 
conditions, based on isotropic lines. 

An isotropic line A is a maximal set of N distinct points in Fq including a — (0, 0) and obeying: 

ct', cr e A ^ (ct', cr) = mod N. (35) 

( The qualification 'mod will frequently be left implicit). It is a fact that each point cr G Fq belongs to at least 
one isotropic line. We mention here some useful properties of such lines, and give further details in Section VF The 
maximality condition allows us to say : 

cr e Fo, (cr, a') = for all cr' e A ^ e A. (36) 

This in turn leads to 

a eX^[N -a]e X, (37) 
where [N — a] is defined in ^ . We also have closure under the group composition law ([51) , © in Fq : 
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cr', cr e A ^ [±2cr], [±3cr], ■ ■ ■ , [a + a'] e X. (38) 

In fact the points {a} of A form an (abclian) subgroup of Fq, of order N, with group composition being (component- 
wise) addition mod N as in eqs. (O, ([6]); the content is the same as in eq. ((38|) : 



cr'.aeX^ g{a')g{(T) = g{[a' + cr]), [cr' + a] G A (39) 

From this group structure we see that if any a, [a' + ct] G A are given, then c' G A is uniquely determined. In case 
g{a) for cr G A is an element of order N, A itself is a cycle generated by a and consisting of the N distinct points 
{(0, 0), cr, [2a], [Sd], • • • , [(iV— 1)(t]}. (However a general A need not be of this form). Examples of such cr are a — (l,p) 
and a — {q, 1). 

For any cr' G Fq, we get a (one-dimensional) unitary irreducible representation (UIR) of A by 

cr G A ^ tJ<'^''"'> (40) 

If cr' G A, this is the trivial UIR . If cr' ^ A,(cr, cr') 7^ mod N for some cr G A, so this is a nontrivial UIR. Hence from 
orthogonality of inequivalent UIR's we obtain: 

uj^"^"'^ = TV if cr' G A, if ct' ^ A (41) 

From the relation ([TO)) for D(cr)'s given earlier it follows that the operators -D(cr), cr G A form a mutually commuting 
set, but they may not form a representation of A. 

Isotropic line marginals condition 

Given A, define 

^YS{^)D{<^) (42) 



N 

ctGA 



Clearly = Px,Tt{P\) = I. Now develop P^: 

= ]^ E S{a)S{a')D{a)D{a') 



(T.cr' 



^E{ E S{a)Sia')r(-'-h{a',a)}D{a") (43) 



7V2 

[ct'+(t]=(t" 



From the subgroup property of A: [cr' + cr] goes over all of A; for given cr, cr" G A, cr' G A is unique. The factor 
^{<j ,a) _ _s^^^ gQ jj-^ ^YiQ last expression, for each cr", {• • • } has exactly iV terms, each ±1. So, since Tr(PA) = 1, 

P| = P\ -4=> P;^ is a rank one projection operator 



^ ^4 E ^('x)^('^')T<^'^^>6(rT',a) = 5(a") 



N 

[<t'+<j]=<t" 



S{cT)S{a') = r<'^'''^>e(CT', cr)S'([CT + ct']), Vcr', cr G A (44) 

If this is obeyed, it implies that {S{a)D{a), cr G A} give a true A^-dimensional UR of A. We have to examine : Can 
these conditions be imposed consistently for all isotropic lines A? If yes, to what extent are the 5'(cr) then determined? 
These are the questions we examine next. 
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V. THE N ODD CASE 

If N is odd so is N'^ and hence from the group structure we have 

cr e Fq or A 3 unique a' £ Fq or A such that a = [2a'] 

i.e., any a = {q,p) = i[2q'], [2p']), unique q',p' G Zat. 

g or p even, < N - I : q' = q/2 < {N - I) /2, p' =p/2<{N - 1) /2 

g or p odd, < A^- 2 : q' = {q + N)/2 > {N + l)/2, p' = {p + N)/2 > (N + l)/2 (45) 
In the relation (|44p . setting a ^ a' to get 

5([2a']) = e(a',a') (46) 

and looking at q/p even/odd we find that S{q,p) are unambiguously given by : 

S{q,p) - (-1)'^' (47) 

Does this obey the condition ([44]) for all A? We find that this is indeed so. The expression (— !)'''' was found by taking 
(7 = tr' in ()44p . Now we put this into that equation with a and cr' independent and ask if it is obeyed, i.e. whether or 
not 

e((7', cr) = Si<7')S{a)Si[a' + a])T<''' -"^ (48) 
holds for all a',(T E \. A key observation here is that since N is odd 

{cr', cr) = mod N (cr', cr) mN 

So the question now is whether 

g^^' i^-j _ (^_2^')g'p'+gp+[g+9']b+p']+9'p-p'9 ^qq'j 

We now check the exponent on the RHS in various cases 

q + q' P + p' exponent rhs 

<A^-1 <A^-1 2q'p +1 

<N-1 >N N{q' + q) (51) 

>N <N~1 N{p'+p) (_i)p'+p 

>iV > +N{q' + q + p' +p) (-1)'?'+9+p'+p+w 

So, comparing with e{a',a), we find that ([50)) holds. 

Thus in the N odd case the marginals conditions for all isotropic lines can be satisfied, all the S{a) are determined 
as above. 

This unique solution is related to the Fourier matrix (actually parity matrix) results. The Fourier operator F on 
Ti. has these actions and properties 

F\q) = \q), F\p) = \[-p]),F^F = FF^ =F^ = 1 (52) 
F\) = \[-q]), F'\p)^\[~p]). (53) 



So i^^ = P = parity operator, which is what we need. From the relation (fT9| 

P\P)(l\ = 7?^ E T-'''P'D{a') (54) 



Ar3/2 

o-'ero 



Set p = q and sum to get 

^ = 1^ E { E ^''^'^''^''nr-^'''D{a') (55) 



Ar3/2 
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Next for P = F , by calculating in the \q') basis 



Tr{PD{a)^) = r''^^ ^ ^-^''<5[,+2,'],o (56) 

and in the \p') basis 

Tr{PD{a)^)^r''P ^ t^^^'W'l.o (57) 

which are necessarily equal. Using the first form , for N odd, the Kronecker delta gives 

q = even = 0,2,4, - ■ ■ , N - 3, N - 1 : q' = [N - q/2] (58) 
q = odd = 1,3,5,--- ,N ^2: q' = {N - q)/2 (59) 

_ jq even : r-^Ptj-J'I^-s/^l = 1 

odd : T--9Pa;-P(^-«)/2 = (-l)P 

(-l)'?f (60) 



Tr{PD{a)^) = 



Hence 



Now 



and hence 



^ = ^ E(-ir^(^) (61) 



N 



M^(0,0) = 5(a) = (-1)9P (63) 

Thus for the case of N odd, we see that there is a unique consistent solution for all signs S{a), such that the marginals 
conditions can be satisfied for all A's. The resulting Wigner phase point operators are the same as those known in 
the literature and are characterised by the fact that the phase point operator at the origin is the parity operator as 
in the continuum case. The existence of a unique square root group element g{<^) for each g{a), guaranteed by N 
being odd, is adequate for this purpose. In particular it has not been necessary to survey in any sense the set of all 
isotropic lines A, their orbit structure under SL{2,'Zn) action (see below) etc. 

VI. THE CASE OF iV A PRIME POWER 

Towards handling the case of general N (essentially even N) we may note the following : Any N can be uniquely 
written as the product of powers of (increasing) primes as : 

k 

N = N1N2 Nk = Y[ Nj 

= pj' , Pj = prime : = 2, p2 = 3, ps = 5, ■ ■ ■ , 
Pj — odd j > 2; and rij = or 1 or 2 ••• . (64) 

If ni — 0, N is odd and then previous results of Section V are in hand. We expect something new to arise only when 
rii > 1. 

We consider the case when is a power of a single prime in the rest of this Section, and turn to the general case 
()64p later in Section VII. Simplifying the notation as much as possible for the moment let us write: 

A = p", p prime, n = 0, 1, 2 • • • . (65) 

(Care will be taken to avoid this prime p being confused with the second entry in the pair {q,p) corresponding to a 
general point a QTq — Zn x Zat). 
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A. Isotropic Lines and 51/(2, Zjv) orbits for N = 

For the isotropic lines we have the foUowing results [2^ : 

1. The total number of isotropic lines A is 

AA=(p"+i-l)/(p-l) (66) 

2. The number JV{a) of isotropic lines passing through a point cr G Fq = Z^r x Z^r is computed as follows. Any 
a £ Zjv can be uniquely written as 

a = ao + aip + a2P^ + h a„_ip"~\ 

a, e {0,1,--- j=0,l,--- (67) 

The 'p-valuation of a' is then the smallest j for which aj is nonzero: 

v{a) = p ~ valuation of a 

~ j such that flo = ai = • • • = ftj-i = 0, aj > 1. (68) 

This definition is unambiguous for a > 1, in particular, we have: 

v{l) = v{2) = ■■■ = v{p-l) = 0; 

v{p) = 1, • • • ; vip") = 2 • • • ; • • • ; v{p^-^) = n - l; 

v{N - 1) = v{p" - 1) = as p" - 1 = (p - 1)(1 +p + p^ + ---+ p"-^). (69) 
We supplement these with the convention 

v{Q) = n (70) 
based on — mod A^. For a = {q' ,p') E I^n 'x '^n define the p- valuation by 

via)^{viq'),vip'))< , q',p'eZN. (71) 

Thus for instance; 

«((0, 0)) = n; z;((l,p')) = v{{q', 1)) = i;((iV - l,p')) = viiq\ N - 1)) = 0. (72) 
Then the number of isotropic lines passing through a G Zat x Zat is 

Comparing with ([55]) we see that A/" = A/'((0, 0)): this is consistent with the condition that any isotropic line A 
must contain (0, 0). For a = (1 or A — l,p'), {q' , 1 or A^ — 1) we have v{a) = 0, A/'(ct) = 1, so only one isotropic 
line passes through each of these points. 

We now turn to the group SL{2,Zi\f) 

S'L(2,ZAr) = {A = ^ " \ a,b,c,d£ZN; ad - be = 1 mod N}- (74) 

and its action on Fq and on the family of isotropic lines thereof. For the case at hand viz. N = p^ the order 
|S'L(2,Zjv)| is given by 

\SL{2,'Ln)\^P^''-\p^-1) (75) 
It acts on the points and isotropic lines in Fq as follows: 

AGS'i(2,Zjv) : cr= (g',p') eFo^CT' = (ag' + 6p',cg' + rfp') eFo, 

: A = {(g',p')} = {{aq' + bp',cq' + dp')}. (76) 

From the latter action one finds that |29| 
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1. The Af isotropic lines, divide themselves into 1 + [n/2] orbits under SL{2,1m) action, where [n/2] is the integer 
part of n/2 . They are denoted by Ok{p"), k = 0,1, ■ ■ ■ , [n/2]. For k < n/2 the orbit contains 

AA(Ofe) = b+l)p""''"' (77) 
isotropic lines, while for fc = n/2 in case n is even we have 

mOn/2) = 1. (78) 

One can easily check in both cases that 

[n/2] 

AA(Ofc)=AA (79) 

fc,0,l,--- 



The largest orbit corresponds to fc = and contains {p + l)p" ^ isotropic lines. 

2. Only the largest orbit Oo(p") has the property that it covers all points in Tq = Z^r x Zjv- The {p + l)p"~^ 
isotropic lines in this orbit are all generated by single generators of order N which may be taken to be (l,p') 
for p' e {0, 1, . . . , TV - 1} and (<?', 1) for q' e {0,p, 2p, 3p, • • • , - l)p}. 



B. Isotropic lines in the 2" case 



Now we specialise to the case N — 2", as otherwise N is odd and then the comprehensive results of Section V are 
available. Thus in a sense this is the most important remaining case. Specialising the above statements ( and further 
quoting from ^2^) we now have: 

1. The total number of isotropic lines A is 

TV = 2"+^ - 1 = 2iV - 1 (80) 



2. If fT is of the form {2j, 2k), then from (f75|) . as v{(j) > 1, the number of A's passing through it is 

7V((T (2j,2fc)) > 3. (81) 

If a is of any of the other three forms {2j, 2k + 1), {2j + 1, 2k or 2k + 1), then as w(cr) = the number of A's 
passing through it is 

7V(cr 7^ (2j,2fc)) = 1. (82) 



3. The A's separate into two types 

Type (a): 3N/2 in number, generated by single generators, and comprising a single (the largest) orbit Oo{2"), 
Type (b) N/2 — 1 in number, involving two generators of orders 2^ and 2^* with both r, s nonzero and r + s = n; 
and comprising all the remaining orbits Ofc(2"), fc = 1, 2, • • • , [n/2]. 

4. Every phase point a lies on (one or more) A's of Type (a). The A's of type (b) cover all the even phase points 
(2j, 2fc) only. 

5. The A's of Type (a) separate further into two subtypes: 

Type (al) containing N A's generated by (l,_po) for po = 0, 1, 2, • • • , — 1. 
Type (a2) containing N/2 A's generated by {qo, 1) for go = 0, 2, 4, • • • ,N — 2 

Therefore each a of any of the three types other than {2j, 2k) lies on a unique A of Type (a) according to the 
pattern : 

CT = (2j + 1, 2fc or 2fc + 1) Type (al) 

a = {2j, 2k + 1) Type (a2) (83) 
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C. Marginals property for Isotropic lines in Z2" x Z2" 



The condition that the average of the phase point operators along an isotropic hne A be a one dimensional projector 
is given in l^^ . The other essential conditions on the signs S{a) are the reflection symmetry (I^SI) and the standard 
marginals conditions (PT|) . ([321) • We know from ([57)1 that cr € A implies [TV — a] G A as well. Applying (Hi)) to such 
pairs of points on A's of Type (a), and remembering that any a lies on such a A, we find that the property (p5| follows. 
( For A's of Type (b) this is only partially true as they cover only the even phase points (2j, 2k)). Thus we begin by 
imposing only the requirements ([3T|) . (|32|) . (|44|l on S{a), for all A's of Type (a). 

We see from ([T^ . A'' being even, that 

e{a,a) = l (84) 

Therefore setting a' — <t in (j44|l leads to 

5((2j,2fc)) = l (85) 

This leaves S{{2j + 1, 2fc or 2k + 1)) and S{{2j, 2k + 1)) to be analysed. Each a of the former type is on a unique 
Type (al) A, while each a of the latter type is on a unique Type (a2) A. We apply (UH) in these cases, choosing 
CT = (To = (l,Po) OT (qo, 1) and a' ~ [2jo'o] or [2^0-0] respectively, thus reaching all points a other than (2j, 2fc), and 
relating S at such points to 5(cto): 

S(i23 + 1, m + l)po])) = (-l)(2^Po-[2Jf«l)/^ . 5((1,po)) X 

f lifpo + [2jpo]<iV-l, . .gg. 
|-lifpo + [2m]>^ ' ^ ' 

5(([(2fc + l)go],2fc + l)) - (-l)(2fe'?o-[2fc9o])/iV.^((^^^l)) X 

f lif qo + [2fcgo] < A-1, 

|-1 if go + [2fcgo] ^ 

In the former relation, the choice of 2k or 2fc + 1 determines po uniquely ; in the latter, that of 2j determines 
uniquely. For go — 0, eq. (pij) determines 5((0, 1)) = 1; for po — 0, eq. ([31]) determines ^((IjO)) — 1. The remaining 
37V/2- 2 undetermined signs are S'((go,l)), go = 2,4, • • • ,iV- 2 and 5((l,po)), po = 1, 2, • • • , - 1. 



These conditions may equivalently be written as 
5((2j + 1, [(2j + l)po])) = 5((l,po))x 



1 if ((2j + l)po - {{2j + l)po])/A^ is even 
1 if ((2j + l)po - m + l)po])/A is odd 



5(([(2fc+i),o,(2.+i)) ^ 5((,o,i))x • (^9) 

1 - 1 it ((2fe + l)(jo — [(2fc + l)(7oJ)/A' is odd 
Thus, for n = 1, 2, when N = 2 and 4 respectively, the free signs are indicated thus: 

S{q,p) 

1 1 (90) 
1 5(1,1) 

S{q,p) 

11 1 1 

1 S'(l,l) 5(1,2) 5(1,3) (91) 
1 5(2,1) 1 -5(2,1) 
1 5(1,3) -5(1,2) 5(1,1) 

In summary: the marginals conditions can be consistently imposed on all isotropic lines of Type (a) comprising the 
largest orbit but leaving 3 x 2"^^ — 2 of the S{a) unfixed. 

We next see by low dimensional examples that these conditions cannot be consistently extended to include isotropic 
lines of type (b). For n = 1,N — 2, there are no isotropic lines of type (b). For n = 2, N = A, there is one isotropic 
line of type (b), generated by (2,0) and (0,2). Condition ([i^ when applied to this isotropic hne gives 5(2,2) = — 1 
conflicting with 5(2, 2) = 1 obtained from isotropic lines of type (a). For n = 3, A^ = 8, there are three isotropic lines 
of type (b) generated respectively by {(2,0), (0, 4)},{(0, 2), (0, 4)},{(2, 2), (0,4)} . Again, as for n = 2, A^ = 4 one finds 
that the results of for isotropic lines of type (b) conflict with those for isotropic lines of type (a)- one can not 
impose marginals property on all isotropic lines consistently. 



12 



VII. ISOTROPIC LINES AND ORBITS IN THE GENERAL CASE 

Turning now to the case of a general N we note that ring Z^r can be factored as 

Zn ~ ^Ni X ZjVa X • • • X Zn^ (92) 

The exphcit correspondence between elements of Z^r and those of the rings ZjVj is provided by the Chinese remainder 
theorem which tells us that an element q G Zpf can be uniquely decomposed as 

k 

1 = ^1r ^3 ■ t^j (93) 

where qj = [q mod Nj] g ZjVj , I'j = N/Nj and fij denotes the (multiplicative) inverse of i/j in Zat^. Thus each element 
q G Zpf can be uniquely represented as an array 

— ' {qi,q2,- ■ ■ ,qk}, qteZN, (94) 

In particular the elements and 1 are represented by 

0^{0,0,--- ,0}; 1^{1,1,--- ,1} (95) 
Further, this correspondence has the nice property that 

q + q < — > {91 + 91, 92 + • • ■ ,qk + q'k}, qi and q[ e z^^ 

(96) 

qq' < — ' {qiq[,q2q2: ■■■ , qkq'k}, qi and g- e Zjv, 

In view of this and the properties ([95]) and (|97p we have the following results: 

• A point (7 G Ztv X Z^r can be represented as 

cr < > {CTi,f72, • ■ ■ ,crfe}, (Ti e Zjv, X ZjVi (97) 

• The symplectic product (it, cr') vanishes if and only if each of the components {ai,al) vanish. 

• The group SL{2,Zn) also factorises as 

SL{2,Zn) = SL{2,Zn,) X SL{2,Zn.,) x S'i(2,ZjvJ (98) 

This can easily be seen by considering the case N — N1N2 and verifying that any matrix A € SL{2, Zm) ■ 

a b 
c 



A = [ ^ ^ ] ; ab — cd — I; a,b,c,d G Z„; (99) 



can be decomposed as ^ = A1A2 where 



(ai,l) (&i,0)\ , _fa,a2) (0,^2) 



= ^ Uio) (V;,i) ) ^ sLi2,z,,)- A2 = ^ f^;,- j e sm^; doo) 

From these considerations it is evident that the isotropic lines in ct G Zjsi x Zjsi and SL{2, Z^r) action are completely 
determined by those in each of the factors Z^j x Z^r^. . 

VIII. TOMOGRAPHIC RECONSTRUCTION FOR AT = 2" 

From the discussion towards the end of Section VI it is evident that for the case when iV = 2" we can only insist 
on marginals property restricted to the isotropic hues of Type (a) constituting the largest orbit under SL{2,Zm) 
action. For each choice for the free signs, we can associate with each such isotropic line a rank one projector P\. 
Each isotropic line generates A'^ — 1 other lines 'parallel' to it obtained, for instance, by shifting the points on it by an 
amount (0, i) in the case of isotropic lines of Type al and by an amount (z, 0) in the case of isotropic lines of Type a2 
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with i taking values 1, 2, • • • , A'' — 1. Denoting by (A, i) the hnes parallel to the line A, the projectors P\^i associated 
with them are obtained by the conjugate action of the appropriate displacement operators on P\: 



From these N projectors associated with each isotropic line we can construct — 1 traceless operators T(^\,i) = 
P{x,i) ~ In/N . Each line of Type (a) gives us iV — 1 T's and since there are 3-/V/2 lines of Type (a), we have a 
collection of 3N{N— l)/2 traceless hermitian operators. Given a density operator for an A^-state system, the operator 
Pn — In /N belongs to the N'^ — 1 dimensional real Hilbert space oi N x N traceless Hermitian matrices. The question 
concerning the tomographic reconstruction of p then reduces to the question whether or not the collection of the 
T's above spans the — 1 dimensional real vector space of traceless hermitian operators. This can be checked by 
examining the rank of the Gram matrix associated with the T's. For iV = 2, 4 we have explicitly checked that the 
corresponding Gram matrices indeed have ranks 3 and 15 respectively. Thus, it would seem that even with restricted 
marginals property the construction developed here permits a tomographic reconstruction of the state of an TV-level 
system though in a non optimal fashion - we have {N — 1){N — 2)/2 more T's then the N'^ — 1 required. 



We have shown how to set up Wigner distributions for finite even dimensional quantum systems working entirely 
with a, N X N lattice instead of a 2A^ x 2A^ grid as was found necessary in the existing formalisms. The Wigner 
distributions thus obtained are consistent with a restricted marginals property and are characterised by 3.2"^^ — 2. 
undetermined signs where n is the exponent of 2 in the decomposition of N into prime factors. As a result, for 
instance, for A^ = 2, 4, 8, 16 there are 2, 2^, 2^", 2^^ different possible definitions of Wigner distributions. 

As a curiosity, in the spirit of the work in [2^ . we have also examined the dependence of the eigenvalues of the 
phase point operators for A^ = 2, 4, 8, as a function of the signs that remain free. (For this purpose it is sufficient to 
look at the eigenvalues of 14^(0, 0)). We find that: 

For N — 2 there is only one free sign, S'(l, 1), and the spectrum of W{0, 0) is the same for 5(1, 1) = ±1 

For A^ = 4 there are three distinct spectra for 14^(0,0) depending on the values of the four free signs, S{1, 1) = 
a, S{1, 2) = 6, 3(1, 3) ee c, S{2, 1) = d. They are 

• ((1 + V6)/2, (1 - V6)/2, -1/2, 1/2) corresponding to a = c, 5 = -d and a = -c, 6 = d. 

• ((1 + 2V2)/2, -1/2, (1 - \/2)/2, (1 - V2)/2 corresponding to a = c = 1, 6 = d = l,a = c = -1, 6 = d = -1 and 

a = — c = 1, 6 = —d 

• ((1 + \/2)/2, (1 + \/2)/2, (1 - 2V2)/2,-l/2) corresponding toa = c=l,6 = d = -l,a = c = -1, 6 = d = 1 and 

a = — c = —1,5 = —d 

For A^ = 8, one has 4 and A^ = 16, one has 15 distinct spectra. Thus, although the number of different Wigner 
distributions based on choices for the signs for N = 2, 2^, 2^, 2* is 2, 2^, 2^*^, 2^^, those which have distinct spectra are 
only 1, 3, 4, 15 in number. ( It seems that the number of distinct spectra for N = 2^ equals 2^~-^ if n even and 2" — 1 
if n odd) The question as to what bring about this enormous reduction is under investigation. Further, it would be 
interesting to see if the square root idea developed here works in the case when the coordinates take values in a finite 
field Q and to see how it relates to Wigner distributions in the more general setting based on the theory of frames 
employed in the work of Ferris and Emerson [2^ . 
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^(O, i)PxD'f{0, i) if A is of Type al 
D{i,0)PxD^i,0) if A is of Type a2 



(101) 
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